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Abstract. - We compare numerical and experimental results exploring the behaviour of 
liquid drops moving across a surface patterned with hydrophobic and hydrophilic stripes. A 
lattice Boltzmann algorithm is used to solve the hydrodynamic equations of motion of the drops 
allowing us to investigate their behaviour as the stripe widths and the wettability contrast are 
altered. We explain how the motion of the drop is determined by the interplay between the 
driving force and the variation in surface force as the drop moves between regions of different 
contact angle and we find that the shape of the drops can undergo large periodic deviations from 
spherical. When compared, the numerical results agree well with experiments on micron-scale 
drops moving across substrates patterned by microcontact printing. 



Introduction. - The question of how liquid drops wet and move across a solid surface has 
long caught the interests of academic and industrial communities alike, with applications rang- 
ing from microfluidic devices to ink-jet printing and surface coating. Though much progress 
has been made since the first pioneering work by Young and Laplace, many interesting, unan- 
swered questions remain. One which has recently come to the fore because of experimental 
advances allowing the fabrication of surfaces with mesoscopic hydrophobic and hydrophilic 
regions is the behaviour of drops on chemically patterned substrates. Several authors [1-6] 
have shown that the wetting behaviour on these substrates can be very rich, with the drop 
shapes depending sensitively on parameters such as the dimensions and contact angles of the 
patterning. In this letter we build on this work to address the dynamics of drops moving 
across an array of alternating hydrophobic and hydrophilic stripes, focussing on the centre of 
mass motion as well as the morphological transitions induced by the imposed external flow. 
The drop is pushed by a constant gravity-like acceleration as opposed to [7] where a thermal 
gradient is applied to generate the drop motion. 

(*) E-mail: halimOt hphys.ox.ac.uk 

(**) Present address: Laboratoire de Physico-Chimie des Polymeres, Universite de Mons-Hainaut, 20, Place 
du Pare - B-7000 Mons, Belgium 

(***) Present address: Institute of computational science, ETH Zurich, 8092 Zurich, Switzerland 



© EDP Sciences 



2 



EUROPHYSICS LETTERS 



We focus on three cases. First we consider the simplest case of a drop moving across 
a boundary between stripes of different contact angle which are much wider than the drop 
radius. This allows us to demonstrate the different roles of the applied force and the surface 
force which arise from the change in free energy as the drop crosses a stripe boundary. We 
then look at stripes of width of order the drop radius and choose the simulation parameters to 
reproduce as closely as possible the experimental values. We find pleasing agreement between 
simulations and experiments. We next consider a substrate on which there is a much larger 
difference in contact angle between the hydrophobic and hydrophilic regions. The shape of 
the drop shows strong oscillations as it is pushed across this surface, and spreading across the 
surface perpendicular to the flow becomes important. 

Equations of Motion. - The equilibrium properties of the drop are described by a con- 
tinuum free energy 

^= f {Mn)^^{danf)dV^ f Mns)dS. (1) 
Jv ^ Js 

ipb{n) is a bulk free energy term which we take to be [8] 

Mn) = Pcii^n + l)'(^n - 2^n + 3 - 2(3r^) , (2) 

where Un = {n — ric) /uc, r^j = {Tc — T)/Tc and ric^T^Tc andpc are the local density, critical 
density, local temperature, critical temperature and critical pressure of the fluid respectively. 
/3 is a constant typically chosen to be 0.1. This choice of free energy leads to two coexisting 
bulk phases of density nc(l ± V Pt^ ) • The second term in Eq. models the free energy 
associated with any interfaces in the system, is related to the surface tension via 7 = 
{4:^2K,pc{(3rw)^^'^nc) /3 [8]. We take the third term as the simple Cahn expression V^s(n) = 
—(pUs where Ug is the value of the density at the surface [9]. The surface field (j) can be related 
to the contact angle of the drops on the surface by 



(j) = 2l3Tw\/2pct<i sign(- - 6>)y cos -(1 - cos -) , (3) 

where a = cos~^ (sin^ 0) and the function sign returns the sign of its argument. 

The dynamics of the drop is described by the continuity and the Navier- Stokes equa- 
tions 

dtn ^ da{nua) = ^ , (4) 
dt{nua) + df3{nUaUf3) = -df3Paf3 + yd(3[n{d(3Ua + daU(3 + 5a(3d^u^)] + naa , (5) 

where u, P, and a are the local velocity, pressure tensor, kinematic viscosity, and accelera- 
tion respectively. The thermodynamic properties of the drop are input via the pressure tensor 
P which is calculated from the free energy 

Pa[3 = {Pb{n) - f (9c,n)^ - K.nd^^n)daf3 + n{dan){d(3n) , 
Vh{n) = pc(^n + 1)^(3^^ -2vn^l- 2(3r^). 

Eqs. Q and are solved using a free energy lattice Boltzmann algorithm [8,10,11]. 

We consider a simulation box of size {L^^ Ly^ L^). The surface at z = is patterned by 
(relatively) hydrophilic and hydrophobic stripes of contact angles Oi and O2 and widths Si 
and 62 respectively. A drop of initial radius R lies on this surface and is pushed along the 
x-direction by a constant acceleration a^^. No-slip boundary conditions are imposed on the 
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velocity field at z = and z = — I. We record the drop shape and the position of the 
centre of mass of the drop as a function of time. Simulation parameters which are used for all 
the numerical calculations are: hz = 0.004, = 1/8, = 3.5, T = 0.4, Tc = 4/7, and = 0.1, 
while those specific to a particular simulation are given in the corresponding figure's caption. 

Experiments. - The experiments were performed by dispensing drops of ethylene glycol 
(volume ~ 2 mm^, surface tension a = 47.7 mN.m~^, dynamic viscosity rj = 16.1 mPa.s, and 
density 1.15 g.cm~^) on surfaces (length = 5 cm, width = 2 cm) patterned by microcontact 
printing. The surfaces were silicon wafers coated by a bilayer of 2 nm of chromium and a 
10 nm layer of gold. First, a PDMS stamp was moulded on a striped master achieved by 
standard lithographic techniques [12]. Then a Immol solution of octadecanethiol in hexane 
was poured onto the stamp. The stamp was applied to the surface and was followed by a quick 
wash with a 1 mmol solution of 1F,1F,2F,2F perfiuorodecanethiol in hexane. This resulted in 
alternating striped layers of octadecanethiol (0.8 mm) and perfiuorodecanethiol (1.0 mm) on 
which the contact angles of ethylene glycol are 80° and 90° respectively. The rms roughness 
of the surfaces after patterning was ~ 2 — 3 nm as determined by AFM in contact mode. 

A CCD camera (25 frames per second) linked to a computer was used to record the 
evolution of the drop after it had slid 3 cm from the top of the substrate. This ensured that a 
stationary regime had been reached. The images were then analyzed using ImageJ software, 
which allowed the position of the centroid of the drop to be followed as a function of time. 

Results. - First, as a simple test case, we present numerical results for two wide stripes 
{Si/R = 82/ R = 200/25 = 8.0) of different contact angles, Oi = 110^ and O2 = 130^. The 
hydrophilic and hydrophobic stripes are labelled 1 and 2 respectively. Figs.Q^a) andQjb) show 
the drop position as a function of time, and its velocity as a function of position respectively. 
Away from the borders between stripes the drop attains a constant velocity as expected. It 
moves faster in the region of higher contact angle. This is because it is subject to a velocity 
profile which increases with height above the substrate; the higher the contact angle the 
further from the surface the drop centre of mass, and therefore the faster the drop moves. 

As the drop moves from the hydrophilic to the hydrophobic stripe there is a pronounced 
dip in its velocity. This is due to the increase in the surface contribution to the free energy 
which can be estimated as 

dvE-, = (72^^ - 72^'' - 7f'' + 7^") y{x)dx (6) 

where y{x)dx is the area of the base of the drop passing the chemical border as it moves 
through dx and 7*^^ and 7*^^ are the solid-liquid and solid-vapour surface tensions. Using 
Young's equation to write this in terms of contact angles gives a repulsive force 

Fs = - ^ =7^''(cos(^2) -cos(^i))y(x) (7) 

where 7^^ is the liquid-vapour surface tension. The interplay between the repulsive wet- 
ting force and the increased efiiciency of the applied force in moving the drop leads to the 
characteristic dip in fig. ^b). If the total force drops to zero we would expect the drop to 
remain pinned at the boundary. Note also that there is a slight increase in the width of the 
drop along the ^/-direction as it reaches the hydrophobic boundary, as it prefers to spread on 
the hydrophilic substrate rather than move across onto the hydrophobic region of the surface. 
However this is disfavoured by the increase in the free energy of the liquid-gas interface and 
is, for the small difference in surface energies considered here, a minor effect. Similarly the 
competition between the attractive wetting force and the decreased efficiency of the applied 
force account for the peak in velocity as the drop passes back into the hydrophilic region. 
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We next consider the case where the stripe width is smaher than the drop diameter, 
(^i/i? = 33/25 = 1.32 and 82/ R = 27/25 = 1.08, and the wettabihty contrast is not too large, 
61 = 80^ and O2 = 90^. These parameters were chosen to match the experimental results as 
closely as possible. The case of a larger wettability contrast will be discussed later. 

We first describe the numerical results. In fig.[2fa) the drop position is plotted as a function 
of time and fig. [2fb) shows the drop velocity as a function of its position on the substrate. 
We also show, in figs.I2fc-h), snapshots of the drop shape at designated points along its path. 
The first difference, when compared to the wider stripes, is that the drop never reaches a 
constant velocity. However, as expected, the periodicity of the substrate is refiected in the 
variaton of the drop velocity with time. The broad features of the curves match those in 
fig. njb) with the fiat sections removed showing that the dynamics are again controlled by an 
interplay between the applied and wetting forces. The additional structure at the peaks and 
troughs of the velocity versus position graph occurs because the drop is traversing two stripe 
boundaries simultaneously. 

We now compare to the experimental results. The simulation parameters were chosen so 
that the values of the liquid viscosity, the surface tension, the drop radius and the stripe 
widths in the simulations matched the experimental values. We also chose the (static) contact 
angles to take values of 80^ and 90^. 

The drop displacement as a function of time measured in the experiments is shown by the 
data points in fig. [^a). Photographs of the drop as it moves across the substrate are shown 
as insets in figs.E^c-h). There is a pleasing agreement between simulations and experiments. 

However we caution that, as with all mesoscale approaches, it is not possible to match 
all the physical parameters and retain a stable and numerically feasible calculation. The 
interface is too wide in the simulations (relative to the drop radius) and the liquid-gas density 
difference is too small. The contact line velocity u ex ^/(An)^, where ^ is the interface width 
and An is the density difference [8]. The drop therefore moves too quickly in the simulations. 
This problem is accounted for by rescaling the time axis. Another difference between the 
simulations and experiments is that experimental drops are not confined in a channel as 
in the simulations. However, we have chosen L^/R sufficiently large that the drop feels a 
parabolic velocity profile akin to the experimental profile. 

However, given these caveats, the agreement between simulation and experiment gives us 
confidence to explore the case of a larger contact angle difference between stripes which is 
more difficult to achieve experimentally. We take Si/R = 0.8, 62/ R = 2.0, Oi = 60^ and 
O2 = 110^. We plot the drop position as a function of time and the drop velocity as a function 
of position in figs.lHfa) and[2fb) respectively. Several authors [4-6] have shown that on such 
a substrate the drop can take two final shapes. One is the 'butterfiy' configuration where 
the drop spans the space between two hydrophilic stripes (eg fig.EIJd)) and the second is the 
'diamond' shape where the drop just lies on a single hydrophilic stripe (eg fig. IHff)). Strong 
deviations from spherical can now be favourable because the extra liquid-gas interface free 
energy is offset by the considerable advantage of lying on a hydrophilic stripe. 

When a constant body force is applied, the drop moves across the surface and changes 
from a diamond to a butterfly shape and back again (flg. IHl^c-f)). Let us assume the drop is 
initially in a diamond conflguration (flg. O^f)). Due to the Poiseuille flow fleld, the drop is 
pushed forward onto the hydrophobic region. Its velocity decreases because of the dewetting 
force at the hydrophobic stripe. If the external body force is small, the drop velocity will fall 
to zero, the drop will be pinned, and the steady-state drop shape will look similar to flg.E^c). 
For the parameters we consider here, however, the drop is just able to channel to the next 
hydrophilic stripe. The effective capillary force at the hydrophobic- hydrophilic border then 
starts to take charge and the drop accelerates and wets the next hydrophilic stripe. The drop 
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(a) Time (10^ simulation time steps) (b) Drop position (lattice spacings) 

Fig. 1 - Drop dynamics on stripes wide compared to the drop diameter: (a) drop position as a 
function of time, (b) drop velocity as a function of position. Simulation parameters: {Lx, Ly, Lz) = 
(400, 100, 80), Oi = llO'^ (dark grey), O2 = 130^ (light grey), 61/R = 200/25 = 8, 62/R = 200/25 = 8, 
and ax = 10~^. 

now has a butterfly shape (fig. Eld)) and is at the peak of an energy barrier between two 
diamond shapes on successive hydrophilic stripes. It then becomes more advantageous for 
the drop to spread along the new hydrophilic stripe in the ^/-direction than to continue to 
move along the substrate. Hence the diamond config uration is re-formed and the oscillations 
repeat. 

Summary. - We have presented numerical and experimental results following the dynam- 
ics of micron-scale drops moving across a chemically patterned surface. The drop behaviour 
is determined by the balance between the driving force and the variation in the surface force 
as it moves between regions of different wettability. As the wettability contrast increases there 
are large deviations of the drop shape from spherical as it moves across the surface. 

We have shown that there is close agreement between results of the lattice Boltzmann 
simulations and experiments. Therefore we hope that the numerical approach will provide 
a useful tool for designing new ways of controlling the motion of liquid drops on patterned 
surfaces for applications in areas such as microfiuidics and ink-jet printing. 

* * * 
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(a) Time (10^ simulation time steps) ^ (b) Drop Position (lattice spacings) 




Fig. 2 - Drop dynamics on stripes smaller than the drop diameter and with a small contact angle 
variation between the stripes: (a) drop position as a function of time. The continuous line shows 
the simulation results and the circles are the experimental results, (b) drop velocity as a function of 
position, (c-h) drop morphology at positions indicated in (b). The large figures are simulations and 
the insets are the experimental results. Simulation parameters: {Lx, Ly,Lz) = (180, 100, 80), Oi = 80^ 
(dark grey), O2 = 90^^ (light grey), 61/R = 33/25 = 1.32, 62/R = 27/25 = 1.08, and = 5 10"^ 
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Fig. 3 - Drop dynamics where there is a large contact angle variation between the stripes: (a) drop 
position as a function of time, (b) drop velocity as a function of position, (c-f) drop morphology 
at positions indicated in (b). Simulation parameters: (Lx, Ly, Lz) = (140,100,80), Oi = 60° (dark 
gray), O2 = 110° (light gray), 61/R = 20/25 = 0.8, 62/R = 50/25 = 2.0, and = 2.5 10"^ 
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